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Abstract: A spatial six degree freedom mathematical model of 
a three wheeled vehicle (TWV) used in Asian countries has 
been developed using multi body system approach. The model 
includes suspension and tyre compliance. The model consists 
of a single sprung mass (vehicle body) connected to three 
unsprung masses (three wheels) at each corner. The 
suspensions and tires are modeled as springs and linear 
damper elements. The model consists of six degrees of freedom 
because the body has three degrees freedom for bounce, pitch 
and roll motions and each unsprung mass have bounce motion. 
Vertical dynamic response of the TWV has been found when 
the vehicle is moving at 45 kniph on random road surface. 
Damping ratios and natural frequencies are obtained using 
Eigen value analysis. Ride analysis has been carried out in 
the frequency domain by performing the spectrum analysis 
using MATLAB/Simulink. 

Index terms: three wheeled vehicle, random road profile, ride, 
frequency analysis, MATLAB/Simulink. 

I. Introduction 

Three wheeled vehicles are extensively used for public 
transportation for small destinations in India and in many 
other countries of Asia. Three wheeled motor vehicles, 
typically used in India and most of the developing countries 
have their front steering with one wheel similar to those of 
motor cycles and two rear wheels are driving wheels with a 
differential and suspension, which are similar to those of 
automobiles. . The term ride is commonly used in reference to 
tactile and visual vibrations. The vehicle is a dynamic system, 
but only exhibits vibrations in response to excitation inputs. 

The response properties determine the magnitude and 
direction of vibrations imposed on passenger's compartment. 
The understanding of ride involves the study of ride excitation 
sources and basic mechanics of vehicle vibration response. 
Ride comfort problem mainly arises due to surface 
irregularities. A mathematical model of a three-wheeled all- 
terrain vehicle(ATV)-rigid rider system without suspension 
with six degrees of freedom (DOFs) has been developed and 
simulation of ATV passing over three bump profiles, of 
rectangular, parabolic, and sinusoidal shapes, has been 
analyzed by Tan and Huston [l].The finite element stress 
analysis of three wheeler chassis [2] has been obtained under 
critical loads, simulating Indian road conditions by 
considering dead weight of vehicle, passengers, driver and 
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ground excitations. Dynamic Vehicle response of a three 
wheeled motor vehicle in frequency domain has been studied 
using finite element modal. The response has been compared 
between the Indian and International roads and different 
modes which affect the passenger comfort [3]. 

A ride comfort simulation model based on the vibration 
of the two-mass system of vehicle body and wheels has been 
build and simulated for the vibration characteristics of the 
model by using simulation software MATLAB/simulink. The 
vehicle ride comfort is evaluated by comparison of the system 
parameters, such as natural frequency of vehicle body, 
damping ratio [4] . Vibration characteristics of vehicle without 
suspension and with front axle suspension were compared 
using 2DOF twin-shaft vehicle dynamic model and 3DOF twin- 
shaft vehicle dynamic model of front axle suspension vehicle 
using MATLAB/SIMULINK, with the white-noise for random 
excitation [5]. The dynamic response of the suspension of a 
road vehicle has been found using experimental setup fitted 
with dampers provided with strain gauges and simulated the 
behavior of the suspension of motor vehicles under the 
control of vibration using a model that more faithfully 
reproduces the actual behaviour. Simulation results in 
MATLAB /Simulink based on the mathematical model 
developed are compared with the experimental data and find 
a good concordance between experimental data and those 
provided by the mathematical model [6] . 

In this paper the mathematical model of TWV with 
suspension using multi body system approach is presented. 
Eigen value analysis has been carried out to find damped 
frequencies and ride characteristics have been studied in 
frequency domain by obtaining dynamic response under 
random road excitation using MATLAB/Simulink. Simulink 
is a software package for modeling, simulating, and analyzing 
dynamical systems. It supports linear and nonlinear systems, 
modeled in continuous time, sampled time, or a hybrid of the 
two. Systems can also be multirate, i.e., have different parts 
that are sampled or updated at different rates. For modeling, 
Simulink provides a graphical user interface (GUI) for building 
models as block diagrams, using click-and-drag mouse 
operations. This is more advantageous compared to other 
simulation packages that require formulating differential 
equations and difference equations in a language or program. 
Simulink includes a comprehensive block library of sinks, 
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sources, linear and nonlinear components, and connectors. 
We can also customize and create our own blocks. Models 
are hierarchical, so we can build models using both top-down 
and bottom-up approaches. After defining a model, we can 
simulate it, using a choice of integration methods. Simulation 
results can be seen while simulation is running using scope 
and other display blocks. 

In addition, we can change parameters and immediately 
see what happens, for "what if" exploration. The simulation 
results can be put in the MATLAB workspace for post 
processing and visualization. Model analysis tools include 
linearization and trimming tools, which can be accessed from 
the MATLAB command line, plus the many tools in MATLAB 
and its application toolboxes. And because MATLAB and 
Simulink are integrated, we can simulate, analyze, and revise 
our models in either environment at any point. 

II. MODELING 

The configuration of six degree of freedom is described 
in fig. 1 




Fig 1 . Discrete element model of the TWV 

Vehicle attitude and trajectory through the course of 
maneuver are defined with respect to a right-hand orthogonal 
axis system, the inertial frame U-X Y Z that is fixed to the 
earth. The origin of the moving reference frame G-xyz coincides 
with the centre of mass of the vehicle body and travels with 
the vehicle. The body of vehicle has 3DOFi.e. rotation about 
x and y axes (roll, pitch) and linear motion along z direction. 
Three wheels have linear motion in z direction (bounce). Three 
independent Euler angles are used to describe the orientation 
of the body-centered frame G-xyz in relation to the iner inertial 
frame. The transformation matrix is developed on the basis of 
these rotations. The Euler angles involve three successive 
rotations about three axes that are not orthogonal in general. 



The concept of yaw, pitch, and roll angles is used while 
selecting axes for rotations. The Newton's second law which 
states that sum of external forces acting on a body in a given 
direction is equal to the product of its mass and acceleration 
in that direction has been used for analysis of the system. 
Force balance equations have been derived for both 
translational and rotational systems. Euler transformations 
are used to simplify the equations. The final six degrees of 
freedom equations that has been obtained are 

The bounce equation of motion for sprung mass: 

fl^z-r [q +c 2 +^]i+ [cjj -c£fr+ [cj/4 +c,/ 4 -cj 3 ]p 

—CjW rr - c^w ri — CM* + [A^ +A^ +k^z+ [kj 2 -kj^ 



(1) 



The roll equation of motion for sprung mass: 

Uf. -+^ 2 y +[c^i -c^-iV + \fiK -+ c i^y-^ c i^\. 
-Cj/jTi;, +[t J /J + -cJk]p+[Kk -Qi¥ +[ j ^f +^i~}' 

The pitch equation of motion for sprung mass: 
[Cj^" + cJ+ + cjfjp -cj^\, + CjAiV,. + [A;£l' 4 - kj^ 4 y 

[V* + V* - *&> + [V/ + W + %C]p- Kh™„ 

- k 3 l$v rl - k^Wj. = 

The bounce equation of motion of rear right wheel: 

m TT w rr -c 1 z +c 2 l 1 r+cjf -c 2 I 4 p+[c 2 +c i ]w„ 

- k 2 z + kj 4 p ■+ [kj + k 5 ] w rr = 

The bounce equation of motion of rear left wheel: 

m A - c i* + c b^ ~ C J*P + ta + c s\™r; 

- k 3 z + k i I 2 r - k^ A p + [k, + k 6 ]w ri = 

The bounce equation of front wheel: 

BIjWj, - Cj Z ■+ cJip 4- [Cj + Cj }l*y ~ Aj Z + A}/jj? 

+ [ky ■+ Ar 4 ]w> = 



■■(:> 



..(3) 



..(4) 



..(5) 



-<«) 



In the equations above z, r, p are bounce, roll, pitch dis- 
placements of sprung mass respectively and w rf w rr> w are 
bounce displacements of rear left, rear right, front unsprung 
masses respectively. The variables that are having single dot 
'.' in super script are velocity components and variables hav- 
ing double dot are acceleration components of respective 
variables. The remaining variables and their numerical val- 
ues that have been taken for simulation latter on in this paper 
are in table 1 . 
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Table I. Numeericle values for simulation 



Parameter 


Variable 


Value 


Units 


Sprung mass 


m. 


492.3 


Kg 


Front unsprung mass 


m.' 


B.5 


Kg 


Rear right unsprung mas s 


1M.T 


5 


Kg 


Rear left unsprunE mass 


m.: 


5 


Kg 


Front damping co ef ficient 


C: 


3500 


N. s,'m 


Rear right damping c o eft 


G 


1107.5 


N. 5 /m 


Rear left damping c o ett 


C- 


2207.5 


N.s/m 


Rear right spring c o eft 


£ 


50400 


N/m 


Rear lett spring co eft 


& 


45400 


N/m 


Front spring coeff 


E 


32700 


Mm 


Front tire damping 


G 


557 


N.s/m 


Rear right tire damping 


C: 


43« 


N. s,'m 


Rear left tire damping 


Ct 


436" 


N. 5,'m 


Front tire stiffness 


E, 


23E2(50 


N/m 


Rear right tire stitthes s 


$ 


250450 


N/m 


Rear left tire stiffness 


s 


250490 


N/m 


Rollins moment of inertia 


4. 


IB2.2 


Kg.nT 


Pitching moment o f inertia 


4. 


ro 


Kg.rn" 


Right wh eel track width from 
central pLane 


I: 


0.575 


m 


Left wheel track width from 
central plane 


h 


0.575 


m 


Distance between front tire 
and C.Ci 


h 


1.49b" 


m 


Distance b etween rear tire and 
C.G 


L 


0.504 


m 


Distanc e b etween C. G- and roll 
axis 


h 


0.396 


m 



IE. RANDOM ROAD EXCITATION 

Road roughness is described by the elevation profile along 
the wheel tracks over which vehicle passes. One of the most 
useful representations is the power spectral density (PSD). A 
plot of the amplitude versus spatial frequency is the PSD. 
The relationship between the power spectral density and 
spatial frequency can be approximated by 



s f (n) = c o- 



m 



Where S o (0) is the power spectral density function, C s 
and N are constants. O is the spatial frequency [4]. 

For vehicle vibration analysis, it is more convenient to 
express the power spectral density of surface profiles in terms 
of the temporal frequency in Hz rather than in terms of the 
spatial frequency, since vehicle vibration is a function of 
time. The transformation of the power spectral density of the 
surface profile expressed in terms of the spatial frequency 
S o (0) to that in terms of the temporal frequency S (f) is 
through the speed of the vehicle v. 



■(f) 



V Q > = - 
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IV. SIMULATION 

Simulation has been carried out in MATLAB/Simulink 
when the vehicle is moving at 45 Kmph on high way with 
gravel road surface (C =4.9xl0 6 andN=2.1). In majority cases 
dynamic systems (that are continuous in time) will be 
described by differential equations. Thus in simulink we 
describe the system with a block diagram and simulate the 
reaction of the system to an input signal. 

Road profile has been described as a PSD function 
(frequency domain) and the simulink has more emphasis on 
dynamic systems (time domain). Hence the road has to be 
generated as a random signal in time domain. The longitudinal 
positions of the signal can be related to a vertical 
displacement of the random road profile at that particular 
point. Thus the simulation of vehicle response on simulink in 
involves generation of random road signal, constructing the 
model and running the simulation for desired time and finding 
the response in frequency domain. 

Random road profile has been generated in time domain 
in MATLAB workspace using sinusoidal approximation 
method [8, 9] in which a single track can be approximated by 
a superposition of N sine waves. 



d (f) = y] A n sin{ nGit-jt - 






&l 



m 



Where tne Fundamental temporal trequency 
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The signal that has been generated from equation (9) has 
been given as input signal to the model. The simulink model 
has been constructed using block libraries. Each mode has 
been modeled separately in different sub systems and they 
are connected using bus creators and bus selectors. The 
parent diagram simulink modal is as shown in fig 2. The 
random road disturbance has been given as input to the 
system at three wheels. The accelerations of different modes 
are obtained in time domain. 

The responses obtained in time domain are exported to 
MATLAB workspace and they are transformed into 
frequency domain using fast Fourier transforms. 






: -F] z J] E t 



i 



-&. 






Fig 2. Parent diagram of simulink model 
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Fig 3. Simulink model for sprung mass bounce mode 




Fig 4. Simulink model for rear left unsprung mass system 

RESULTS AND DISCUSSIONS 

In the present work ride characteristics has been studied 
for the TWV in frequency domain when the vehicle is moving 
at 45 kmph on with gravel road profile. Natural frequencies 
have been found using Eigen value analysis. Out of twelve 
Eigen values that have been obtained from Eigen value 
analysis, two of them are found to be real and negative. . The 
remaining ten are complex occurring in conjugate pairs 
representing five oscillating modes. Out of five oscillatory 
modes of vibration, three of them are identified as roll, pitch 
and bounce modes for sprung mass system and remaining 
two are for unsprung mass systems. The damped frequencies 
and corresponding damping ratios for different modes that 
have been obtained from Eigen value analysis are in table 2. 

Table II. Damped frequencies of significant modes 



S.NO 


MODE DES CAP- 
TION 


DAMPED FRE- 
QUENCY (HZ) 


DAMPING 
RATIO 


1 


Sprung mas s b ounce 


2.47 


0.35 


2 


Sprung mass roll 


1.73 


0.24 


3 


Sprung m a; s pitch 


5.94 


0.64 


4 


Rear right unsprung 
mas 5 bounce 


17.21 


0.B04 


5 


Rear left unsprung 
mas s b ounce 


17 21 


0.B04 


6 


Front unsprung mass 
bounce 


is.:: 


0.829 



The response of the vehicle i.e. power spectral densities 
of acceleration of heave, roll, pitch of sprung mass are as 
shown in fallowing figures. 
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Fig 5. PSD of Acceleration of sprung mass bounce 
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Fig 6. PSD of acceleration of sprung mass roll 
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Fig 7. PSD of acceleration of sprung mass pitch 
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Fig 8. PSD of acceleration of rear right unsprung mass 
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Fig 9. PSD of acceleration of rear left unsprung mass 
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Fig 10. PSD of acceleration of front unsprung mass 

From the results of sprung mass we can easily observe 
enhanced response at certain natural frequencies. Roll is 
having less impact on vertical vibrations compared to that of 
pitch. The acceleration of pitch has its maximum at 6 Hz and 
we can see increase in sprung mass bounce at the same 
frequency. Pitch mode is a more annoying when compared 
with other modes because its effect directly transferred to 
the passenger in the vehicle and they feel discomfort. PSD of 
acceleration of bounce of sprung mass has a maximum value 
at 2.8Hz. There a considerable increase in the acceleration in 
the range of 4Hz-8Hz which is more annoying because 
passengers in the vehicle feel discomfort because in this 
frequency range the resonance of the human stomach occurs. 
The unsprung mass natural frequencies are about 8 times to 
that of sprung mass bounce frequency. The responses of 
unsprung masses are changing in entire frequency range and 
wheels are less damped compared to the sprung mass. Wheel 
natural frequencies are high they have less impact on human 
comfort. Ride assessment using ISO tolerance limits are as 
shown in figll. 
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CONCLUSION 

A spatial six degree freedom of TWV has been developed. 
Vertical dynamic response of vehicle has been obtained when 
the vehicle is moving at 45 kmph on random surface using 
MATLAB/Simulink. Ride analysis has been carried out in 
frequency domain. Damped frequencies and damping ratios 
has been found using Eigen value analysis. 
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Fig 11. Comparison of vehicle vertical vibration with ISO tolerance 

limits 
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